clear

load('kasagi.mat')
load('periodic.mat')
load('recycling.mat')
load('jicf.mat')

x_range=[0 150];
linear_log='linear';
%linear_log='log';

figure1=figure;
plot(kasagi.y,kasagi.u_mean,'LineWidth',2);
hold('all');
plot(periodic.yp,periodic.um,'.','MarkerSize',15);
plot(recycling.yp,recycling.um,'+','MarkerSize',10);
plot(jicf.yp,jicf.um,'x','MarkerSize',10);
figure_property( figure1, 'u\_mean vs y+', 'y+', 'u\_mean', x_range, linear_log );
legend('Kasagi DNS','Periodic','Recycling','Recycling+jet','Location','SouthEast');

if (true)
figure2=figure;
plot(kasagi.y,kasagi.uu,'LineWidth',2);
hold('all');
plot(periodic.yp,periodic.Ruu,'.','MarkerSize',15);
plot(recycling.yp,recycling.Ruu,'+','MarkerSize',10);
plot(jicf.yp,jicf.Ruu,'x','MarkerSize',10);
figure_property( figure2, 'u\_rms vs y+', 'y+', 'u\_rms', x_range, linear_log );
legend('Kasagi DNS','Periodic','Recycling','Recycling+jet');

figure3=figure;
plot(kasagi.y,kasagi.vv,'LineWidth',2);
hold('all');
plot(periodic.yp,periodic.Rvv,'.','MarkerSize',15);
plot(recycling.yp,recycling.Rvv,'+','MarkerSize',10);
plot(jicf.yp,jicf.Rvv,'x','MarkerSize',10);
figure_property( figure3, 'v\_rms vs y+', 'y+', 'v\_rms', x_range, linear_log );
legend('Kasagi DNS','Periodic','Recycling','Recycling+jet','Location','SouthEast');

figure4=figure;
plot(kasagi.y,kasagi.ww,'LineWidth',2);
hold('all');
plot(periodic.yp,periodic.Rww,'.','MarkerSize',15);
plot(recycling.yp,recycling.Rww,'+','MarkerSize',10);
plot(jicf.yp,jicf.Rww,'x','MarkerSize',10);
figure_property( figure4, 'w\_rms vs y+', 'y+', 'w\_rms', x_range, linear_log );
legend('Kasagi DNS','Periodic','Recycling','Recycling+jet','Location','SouthEast');

figure5=figure;
plot(kasagi.yt,kasagi.t_mean,'LineWidth',2);
hold('all');
plot(recycling.yp,recycling.tm,'.','MarkerSize',15);
figure_property( figure5, 't\_mean vs y+', 'y+', 't\_mean', x_range, linear_log );
legend('Kasagi DNS','Recycling','Location','SouthEast');

figure6=figure;
plot(kasagi.yt,kasagi.t_rms,'LineWidth',2);
hold('all');
plot(recycling.yp,recycling.t_rms,'.','MarkerSize',15);
figure_property( figure6, 't\_rms vs y+', 'y+', 't\_rms', x_range, linear_log );
legend('Kasagi DNS','Recycling','Location','SouthEast');

figure7=figure;
plot(kasagi.yt,kasagi.vt,'LineWidth',2);
hold('all');
plot(recycling.yp,-recycling.vtm,'.','MarkerSize',15);
figure_property( figure7, '< -vt > vs y+', 'y+', '< -vt >', x_range, linear_log );
legend('Kasagi DNS','Recycling','SouthEast');

end


